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Abstract 

The piston problem is investigated in the case where the length of the cylinder is infinite (on 
both sides) and the ratio m/M is a very small parameter, where m is the mass of one particle of 
the gaz and M is the mass of the piston. Introducing initial conditions such that the stochastic 
motion of the piston remains in the average at the origin (no drift), it is shown that the time 
evolution of the fluids, analytically derived from Liouville equation, agrees with the Second Law 
of thermodynamics. We thus have a non equilibrium microscopical model whose evolution can be 
explicitly shown to obey the two laws of thermodynamics. 

Keywords: Noncquilibrium, thermodynamics, entropy, entropy production, piston, similarity flow. 

1 Introduction 

Recently the well-known "adiabatic piston problem" has attracted a lot of attention [l]-[25]. Most 
investigations have however concentrated on the motion of the piston and very few studies have been 
made concerning the time evolution of the fluids on the two sides of the piston. In our previous works 
[5, 6] it was shown that the piston in a cylinder of finite length evolves with two different time scales. 
In a first time scale the motion of the piston can be considered as "deterministic" and "adiabatic" , 
and the system evolves to a state of mechanical equilibrium where the pressures on both sides are 
approximately equal, but the temperatures different. Then on a much larger time scale (if m <C M), 
the piston evolves "stochastically" and with "heat transfer" to a state of thermal equilibrium, where 
the temperatures (and the pressures) on both sides of the piston are equal. It was shown in [5 a] that 
in the first time scale the relaxational motion of the piston can be either weakly damped or strongly 
damped depending on whether the parameter R = M gas /M pis ton is small (R < 1) or large (R > 1) 
where M gas is the total mass of the gas. 

A microscopical analysis of this model was given in [5 b]. We were able to derive equations 
describing the time evolution based on three "reasonable" assumptions. It was also shown that the 
results obtained from these equations were qualitatively in good agreement with those observed on a 
very large number of numerical simulations, except for the observed damping coefficient which appears 
a lot smaller than the predicted value. It was realized that the underlying damping mechanism must 
be related to the propagation of sound waves bouncing back and forth between the boundaries of the 
cylinder and the piston. Unfortunaly the effect of such a wave propagation could not be taken into 
account in our previous analysis, because of the "average assumption" we had introduced (i.e. the 
values of density, pressure, temperature on the right/left surface of the piston can be replaced by the 
average of these quantities in the right/left compartment). Therefore to understand how the damping 
effect can be generated by purely elastic collisions on the piston surface and cylinder walls, we are 
forced to investigate the propagation of wave through the ideal and out-of-cquilibrium fluids. 

To simplify the following investigation on the wave propagation in the fluids, we shall consider 
the case where the length of the cylinder is infinite on both sides and all velocities (of the particles 
and the piston) are parallel to the axis of the cylinder. Initially the fluids on the left and on the 
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right of the piston are each in thermal equilibrium with densities nj, temperatures T ± (T + > T ~) 
and pressures pf , where the index +/— refer to the right/left of the piston. To simplify the problem 
further we have choosen the pressures p$ (given TJ^) in such a manner that the piston, which evolves 
stochastically under the elastic collisions with the particles, remains on the average at the origin (for 
a time sufficiently short to neglect the Brownian deviation - a second-order contribution to the piston 
motion). With these initial conditions no work will be done by the right side on the left side and there 
will be only heat transfer (from the right to the left since Tq~ > Tq). Following the analysis in [5 b], 
which is a singular perturbative approach in powers of the small parameter m/M, we shall derive, for 
these initial conditions, the distribution function for the velocity of the piston, in the stationary state 
(of the piston), to first order in ^JrajM . Using the fact that we have only elastic collisions in our 
model, we shall then derive the distribution function p(x, v; t) for the fluids. Finally we shall obtain 
from the well-known relations of kinetic theory the density fields n(x, t), the velocity fields w(x, t), the 
stress (tensor) r(x,t) and the "heat" flux jQ(x,t). In other words we will obtain the time evolution 
of our microscopical model, considered as a fluid out of equilibrium. 

The next problem we shall discuss in the following is whether the evolution thus obtained obeys 
the laws of non-equilibrium thermodynamics. In particular the question will be whether it is possible 
to find a microscopically rooted entropy function s from which, following the laws of thermodynamics, 
we recover the evolution of the fluids. 

It will be shown that is possible to consider our fluid as a two-component simple fluid described by 
a fondamental equation s = s(u, m, n 2 ) from which we can recover our fluid equations, together with 
the well-known phenomcnological equations of thermodynamics, and an entropy production which is 
strictly positive. Moreover this entropy function will coincide with the Boltzmann entropy at those 
points where the fluid is in a stationary state (i.e. at ±00 and on the right/left surface of the piston). 

In Section 2 we shall recall the thermodynamic equations of the I-componcnt and the 2-component 
simple fluids. We then recall in Section 3 the microscopical definitions introduced in kinetic theory 
for the thermodynamical variables. The piston problem and the distribution function for the velocity 
of the piston in the stationary state are presented in Section 4. We then discuss the properties of the 
fluids at the surface of the piston in Section 5, and the general evolution of the fields in Section 6. 
The problem of entropy and entropy production is discussed in Section 7-8 and numerical graphs of 
the different fields are presented for a specified set of initial conditions. Finally, general conclusions 
are presented in the last section. 



2 Thermodynamics of 1-dimensional simple fluids 

It is well known that the conservation of the number of particles, the linear momentum, and the 
energy, i.e. the First Law of thermodynamics, give for a single component fluid in I dimension the 
following equations for the density of particles n(x, t), the fluid velocity w(x, t) and the energy density 
e(x, t): 

d t n + d x (nw)=0 (1) 
mn[dtw + wd x w] — d x r = (2) 
d t e + d x (ew + j e ) = (3) 

where m is the mass of one particle, r is the stress (tensor) and j e is the energy current, adding to 
the "convective" contribution ew. We have also assumed that the fluid is not submitted to external 
forces. Furthermore, the Second Law of thermodynamics yields for the entropy density s(x,t): 

d t s + d x (sw + j s ) = i 

(4) 

i = i(x, t) > for all (x, t) 
where js{x, t) is the entropy current and i(x,t) the "internal entropy production" or "irreversibility". 
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A simple fluid is denned by the condition that there exists a function u — u(s,n), independent of 
x and t, such that: 

e = - mnw 2 + u(s, n) (5) 

and in this case, one introduces the temperature T(x,t) and the pressure p(x, t) by the thermody- 
namical definitions: 

_ du du m 

T=— , fj, = g^, p = Ts + nn-u (6) 

It follows from Eqs. (1-6) that: 

j e = -tw + jq (7) 
3Q = Tjs (8) 

* = Jgdx (Jj^J + ^(r+p)d x w (9) 

The description of a two-component fluid involves an additional equation, for the density m(x,t) 
of the first component and the associated current j\{x,t): 

d t m + d x {n lW + =0 (10) 

Assuming that the particles of the two components have all the same mass, the simple fluid is now 
defined by the condition that there exists a function u — u(s, n\, n-i) such that: 

e = - mnw 2 + u(s, m, n 2 ), with n = n\ + n 2 (11) 

Again one introduces the temperature T(x, t) and the pressure p(x, t) by the thermodynamical defi- 
nitions: 

—. du du 

T=—, = — , p = Ts + mm + fi 2 n 2 - u (12) 



In that case we obtain 



j e = -TW + j Q (13) 

jq = Tj s + ji(m -m) (14) 

i = J Q d x +hr+ p)d x w + ji d x (^-A) (15) 



Let us stress that the "heat current" is not a uniquely defined concept in the case of a two-component 
fluid. On can either introduce the heat current by the condition that the sum of the "heat" current 
plus the "work" current —tw yields the "energy" current j e , so that j'q defined above by Eq. (13) is 
interpreted as the heat current. Alternatively one can define a heat current ju by the condition that 
Eq. (8), written above for a one-component fluid, remains the phenomenological relation between the 
"heat" and the "entropy" currents, i.e. 

3h = Tj s = Jq - ji{m - M2) (16) 

(17) 

With this definition the entropy production is expressed by 

i = 3h9x (Jj^J +j;(r+ p)d x w + ^ ji d x (u 2 - m) (18) 
which is the usual form in thermodynamics. 
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The irreversibility i(x,t) is thus a sum of products ["current" times "force"] and usually one 
introduces at this point phenomenological relations between forces and currents to ensure that the 
irreversibility is non-negative for all (x, t). Indeed, whereas the First Law Eqs. (1-3) can be understood 
from any microscopic viewpoint and has the status of an exact and inviolable necessity, the Second Law 
Eq. (4), which is of fundamental importance for non-equilibrium macroscopical physics, is a statement 
concerning concepts, "entropy" and "entropy production", for which there seems to be no agreement 
on a microscopical definition, except for equilibrium states, and whose necessity for non-equilibrium 
situations is even subject to controversy [26, 28]. 

In the following we shall consider this issue for the piston problem described in the introduction. 
For this problem we are led to consider the special case where 



1 , m T 
u=- nk B T = -- 



It then follows that: 



s(u, ni,n 2 ) — nks 
where g = #(711,712) is so far undetermined and: 



1 / , 47ra\ 
- 1 + In 3 

2 \ mn 6 J 



+ g(ni,n 2 ) 



dg 



dg 



p = nksT 1 — m— n 2 — — 



dni 



' dn-. 



In that case the viscous-stress (tensor) r/ r is given by: 

f dg dg\ 

Tf r =T+p= -nk B T[ m- hn 2 - — 

\ am ou2 J 



(19) 
(20) 

(21) 
(22) 



Note that if the two fluids are identical and in equilibrium, then using the definition of Boltzmann 
entropy (Section 3) for s(x,t), we have g(ni,n 2 ) = 0. 



3 Kinetic theory 

To obtain a microscopical definition of the entropy density s(x, t), we consider the kinetic framework 
underlying the above fluid description. We consider a system of point particles, with mass m, in a 
semi-infinite cylinder x < 0. We assume that all velocities are parallel to the axis of the cylinder (i.e. 
the x-axis) and that the particles interact through pure elastic collisions only. In kinetic theory, this 
"one-dimensional" system is described by a distribution function p(x, v, t) solution of the continuity 
equation: 

(d t + vd x )p(x, v,t) = for x < (23) 

together with the boundary condition at x = 0: 



dv p(x, v, t)v = (24) 



It is well known that one recovers the conservation equations for the one-component fluid Eqs. (1-3) 
by means of the following definitions: 

/ + OC 
dv p(x,v,t) (25) 
-00 

1 f +oc 

w(x,t) = — -. 7 / dv p(x,v,t)v (26) 

n{x,t) ./-00 

/ + OC 
dv p(x, v, t)(v — w(x, t)) 2 (27) 
-00 
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/+oo 
dv p{x,v,t){v-w(x,t)) 3 (28) 

e = i nmw 2 + u (29) 
u = -\t (30) 

je = -TW+j Q (31) 

We would like to complete these equations with a continuity equation for the entropy. In other words 
we want to define the fields entropy density, current of entropy, production of entropy, as well as the 
temperature and pressure fields. 

One can introduce for example the Boltzmann entropy: 

/+oo 
dv p(x,v,t)lnp(x,v,t) (32) 
-oo 

which obeys the equation: 

d t s B + d x {sBW + j B ) =0 (33) 

with: 

/+oo 
dv p(x, v, t) [v — w(x, t)] In p(x, v, t) (34) 
-oo 

However we do not have a standard definition of the "temperature" . If we identify j B with the current 
of ss, then to obtain the fluid equations (4) and (8), we are led to define: 

T B (x,t) = jQ [ X ^l and i B (x,t)=0 (35) 

jB{X,t) 

The condition i B = is understood from kinetic theory where the internal entropy production is 
induced by the collision term which is zero in our case. The problem with the above definition of 
temperature is that we do not know whether thermodynamic relations between n, u, s Bl T B such as 
Eq. (6) will be satisfied. 

If on the other hand we adopt the usual definition of temperature as related to the thermal energy, 
i.e. for non interacting particles in 1 dimension: 

u = i nk B T (36) 
3Q = Tjs (37) 



and the relation: 
then we have: 

with i B (x,t) given by: 



i B =d x [§-j B ) (39) 



but we cannot conclude at this point that i B (x, t) will be non-negative. Furthermore we do not have 
in general the thermodynamic relation l/T = ds B /du. 

We could also consider that the fluid defined by p(x, v,t) is a simple 1-component fluid with 
temperature given by the thermal temperature Eq. (36). In that case we have for the thermodynamic 
entropy, pressure, and viscous stress: 



s(x, t) = nki 



1 / , 4ttu , 
- 1 + In + g(n) 

2 V mr A ' 



(40) 
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p = nkpT — n 2 kBTg' (n), hence r/ r = r + p = —n 2 kBTg' (n) (41) 

We would have to find at this point an appropriate function g(n) such that ti(x,t), Eq. (9), is non 
negative. The simplest choice would be g = in order that s coincide with ss at equilibrium. 

In the example presented in the introduction (piston problem, see also Section 4) the entropy at 
point (x,t) must depend on 3 parameters (no,To,Tp) where no, To are the density and temperature 
at x = — oo and Tp the piston temperature. It is thus not possible to consider the fluid as a simple 
one-component fluid described by 2 state variables, i.e. s — s(u,n). However, we shall see that it is 
possible to consider the system as a 2-component fluid with an entropy function: 

s(x, t) = s[u(x, t),ni(x, t), n 2 (x, t)] (42) 

which coincide with the Boltzmann entropy at x = — oo and x — — 0, where the fluid is in a stationary 
state. In that case we are led to introduce a temperature T s by the thermodynamic definition: 

_L = |_ sKni , n2 ) (43) 

a current j\ satisfying the continuity equation Eq. (10) and an entropy flux now given by Eq. (14). 
Again with this choice we have to verify whether i(x,t) given by Eq. (15) is actually non negative. 

The purpose of the following sections is to explore what these different definitions imply for the 
piston problem, and whether we can arrive for this model at a description consistent with the ther- 
modynamics of fluids, i.e. the two laws of thermodynamics. 



4 Piston problem 

We consider an infinite cylinder of area A containing two gases separated by a movable piston. The 
gases arc made of point particles with mass m, while the piston is a rigid solid with mass M 3> to and 
no internal degrees of freedom. The particles and the piston interact through purely elastic collisions 
only. We assume that all velocities are parallel to the x-axis (the axis of the cylinder) which reduces 
the system to one dimension. Initially the piston is at rest and the two gases are homogeneous and 
in thermal equilibrium described by Maxwellian distributions of velocities with temperatures T^, 



hereafter denoted 



y T ±, 



and densities respectively on the left (— ) and on the right (+) of the 



piston. We shall take initial conditions n^, such that the piston evolves toward an equilibrium 
state with (V) = 0, and for clarity we choose > T ~ . Using the results of [5 b], this requires that 
the initial pressures p^ — n^ksT^ satisfy the condition: 



VZ -Po = ^(Po +Po) 




'my 



> 



(44) 



Using a singular perturbativc approach to first order in m/M, it was shown that the piston will reach 
an equilibrium state $(V), with (V) = and temperature Tp where 



k B T P = M(V 2 ) = kp \J TqTq 



(45) 



in a time t\ which is proportional to M/A and can be made arbitrary small. In this equilibrium state 
of the piston, the heat passing through the piston per unit time is given by [5 b]: 



p(+)-(-) 
Q 



Am 8k B . + , _, 
2M V^ (??0+Po) 




+ o 



/my 

\m) 



> 



(46) 
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In the present paper, we take this equilibrium state of the piston and X = as the initial conditions 
for the piston and Maxwellian distribution with for the gases. The piston will henceforth diffuse 
around X = 0, but this diffusive motion is very slow compared to the relaxation phenomena inside 
the gases (damping oscillations) and it can be neglected on the time scale here considered. 

In [3], the velocity distribution function of the piston &(V) in the stationary state was obtained 
to first order in y/m/M for the case p^ = Va ■ The same analysis was conducted for the case p^ ^ p^ 
in [7]. For our considerations where p^ and p$ satisfy the condition Eq. (44), one obtains: 



with 



0P = 



-PpV 2 



M 



jl + iv^VM ai ^faV -iy^VM ai (V^) 3 + o( 



2k B \/T+T„ 2ksTp 



1 o 



1/4 



(47) 



(48) 



Assuming that we can neglect the recollisions of particles on the piston, or better, introducing a 
cutoff in the velocity distribution such as was done in [9]-[ll], [14]-[16] to ensure that this condition is 
satisfied, one can thus compute the velocity distribution of the particles at the surface of the piston. 
The distribution of velocities of the particles before (P(v)) and after (P(v)) a collision on the piston 
are given respectively by: 



and 



P~(t;)~ f dV<p T -{v)$(V) 

J — CO 

on the left side and 

P+(u)~ / dVip T+ (v')$(V) and 

J V 

on the right side where: 

v' = -(1 - a)v + (2- a)V 
V = (1 — a)V + av 



P~(v) 



P+(v) 



dVip T - (v')^(V') 



dVip T+ (v)$(V) 



(49) 



(50) 



2m 



a = 



M + m 



< 1 



(51) 



Introducing the boundary condition at the surface of the piston 



/ 

J — ! 



dvp ± (0,v,t)v = 



(52) 



and using the fact that 3>(V) is peaked in V = 0, of width (fcsTp/M) 1 / 2 far smaller than typical v, 
we obtain the distribution function of the particles on the left surface of the piston: 



where 



p-(0,v,t)=6(v)p-(v)+6(-v)p-(v) 
J- 

- e '~ 



P ( v ) = n o \ — 




-0-v 2 



(53) 



(54) 



P («) = 



1 



-f3'v 2 



7 



with 



0- 



P~ = 



2k B T 
1 + 5- 



(55) 



<T =a{2-a) 



-1=0 



m) >0 



Similarly the distribution function on the right surface of the piston is given by: 

p+(0, v, t) = 6{-v)p+{v) + 6{v)p+{v) 
together with Eq. (54) where (3~ , (3~ and S~ are replaced by 

' /?+ = m 



(56) 



2k B T+ 
1 + 6+ 



(57) 



6+ =a(2- 



O 



(-) 

\m) 



< 



Strictly speaking the distributions (53) and (56) are valid only for \v\ ^> (fceTp/M) 1 / 2 , because of 
the error functions which appear in (49) and (50). Indeed, in (49) and (50, we replaced J_ ca dV by 

6[v] f*™ dV , and J^°° dV by 6[—v] f*™ dV , which is valid only if \v\ is far larger than the width 

of $(V), namely \v\ » \J (V 2 ) = (ksTp/M) 1 / 2 . However in the following, we shall be interested only 
in thermodynamic quantities defined by averages with respect to p+(0,v,t) and since M 3> m = 1 
(typically M = 10 3 — 10 6 in our simulations), the corrections due to error functions will be negligible. 



5 Thermodynamical quantities at the surface of the piston 



Given the distributions (53) and (56), the thermodynamical quantities defined in Section 3 are imme- 
diately obtained at the surface of the piston (x = 0). On the left side we have: 



n- m = \n- (l + -^= 



2u-(0, t) = -T~ (0, t) = X - po (l + Vl + tf-) = Po (l 



m / Tp 



M \T< 



O 



/ray 

Vm) 



v ; k B n-(Q,t) 



1 s _ _ 2k B T - 



(58) 

(59) 
(60) 

(61) 



and similarly for the right side, with corresponding quantities , p^ , Tq and 5+, except for j'q where 
an additional sign change is involved: 



(62) 
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As can be seen from (59), the (mechanical) equilibrium condition of the piston, i.e. the equality of 
the forces on both side (since (V) = 0): 



T~(0,t) = r+(0,t) 



for all t 



(63) 



is to first order in m/M identical with the condition Eq. (44) previously obtained at this order (in the 
framework of kinetic theory). Moreover, from the stationarity condition expressed by 



j Q (0,t) = .? + (0,i) =j Q (0,t) for all t 

we recover at lowest order in m/M our previous results Eqs. (45-46), i.e. 

T P = y/T+To + O (hence T ~ < T P < T+) 



1 



my 
MJ 



For the Boltzmann entropy Eqs. (53) and (32) give: 

-u. «(0,f) 









1- to- 




rn ) 



< 



C(6) 



(64) 

(65) 
(66) 



C(5) 



1 



i + VTTs 



ln(l + J)+ln(l + >/m)-In2 = -^* 2 + o[(^) 



(67) 



where the subscript (+/— ) for right/left have been omitted and the currents jg of Boltzmann entropy 
(34) are 



jf(0,t)= ± - ln(l + ^)n±fc B 



'2k B T= 
mil 



< 



(68) 



Finally integrating the continuity equation (33), the Boltzmann entropy Sg on the right/left com- 
partment satisfies equation : 



dt~" y "' "T±(0,i) 
From (66, 68), the "Boltzmann temperature" (35) is: 

S 1 



T B (0,t) = T(0,t) 



ln(l + <5) VTT5 



= T(0,t) ( l + ^ + O 



/my 

\m) 



(69) 



(70) 



i.e. both definition of temperatures coincide at order m/M on the surface of the piston. If we identify 
Tb(0, t) with the thermodynamic temperature then (69) is the well known expression of the Second 
Law and the entropy production I ± (t) on the right/left of the piston is strictly zero. Thus we must 
conclude that i(x,t) = and the increase of the total entropy is entirely due to heat transfer. 

On the other hand taking the "thermal temperature" T(0, t), Eq. (60), as the thermodynamic 
temperature, we obtain for the entropy production: 



dt St{t) T T±(M) 



= ± 



ln(l + S) 



A 



VT+S-1 



A 



T±(0,t) 



24 T±(0,i) 



jq(0,t) + O 



/ m \ 4 

\m) 



(71) 
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Since jq is order of m/M, we thus conclude that i(x,t) is zero to second order in m/M; but jQ(0,t) 
being negative, the definitions ss(x,t) and T(0,t) are not compatible at order 0[(m/M) 3 ] since they 
lead to an entropy production which will be negative for some (x,t). 

Conclusion 

The distribution functions at the surface of the piston p ± (0, v,t), Eqs. (53) and (56), yield ther- 
modynamic properties which are consistent to first order in m/M with previous results and with the 
laws of thermodynamics. In particular, the forces on both sides of the piston are equal, the "heat 
flux" jq is continuous, the different definitions of "temperature" and "entropy" coincide at this order 
on the surface of the piston. Moreover the viscous stress r + p is zero. If wc adopt the Boltzmann's 
definition of entropy, then the total entropy production: 




is entirely due to heat transfer. However with this definition of entropy, we are forced to take Tb as 
definition of temperature, at least if we want to go at order (m/M) 2 . 



6 Propagation through the fluid: similarity flow 

Since the effect of purely elastic collisions between the particles is identical to the non-interating case, 
we consider that the distribution function p(x, v, t) for the gas is the solution of the continuity equation 
(33), given by: 

p(x,v,t) = p~(x,v,t) + p + (x,v,t) (73) 



(74) 



where 



with: 

p~ (x, v, t) = 6(-x) [9(vt - x)p~ (v) + 9(x - vt)p- (v)} 
p+(x,v,t) = 9(x)[6(x-vt)p + (v) + 6(vt-x)p+(v)] 

p(v) = n a tp(v) with tp(v) = \ — e~@ v 

V 7T 

p(v) = n Q y(v) with !p(v) = J- e-^ 2 ( 75 ) 

1 + V 7T 

777 — 

P=— and /3 = / 3( l + ^ 

with 5 given by Eq. (55) and the subscripts (+/— ) corresponding to (x > 0) rcsp. (x < 0) have been 
omitted. We thus have a similarity flow, i.e. for all field f(x,t) we have: 

f(x,t) = m where £=y (76) 

Notice that £ is negative on the left and positive on the right of the piston. 

From the expression p(x,v,t), Eq. (73), we can compute the thermodynamical quantities defined 
in Section 3. We shall detail the computations only for the left side of the piston (£ < and <5~ > 0) 
and shall drop the (— ) subscript. The different fields arc represented on Fig. (1-5) for special values 
of Uq and T^. From the definition of density n(x,i) — n(£): 

H-°° ri 
n(0 = / dvp(v) + / dvp(v) (77) 
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follows that: 

%n(0 = p(0 p(0 (78) 

Defining £ (in fact £ _ but we have omitted the subscript (— ) in this whole section) by p(£) = p(£), 
which gives (note that £ < 0): 



^ = (i + j) Ml+i) hence f - 



m 



/ m \ 



we have 



p(£)<p(£) for £<£<0 



,?(£)> p(£) for 

and n(£) reaches a maximum at £ = £ (see Fig. 1). For the velocity field we have (Fig. 2): 



n(£M£)=n ^^ ( 



-/3« 2 



(79) 



(80) 



(81) 



Therefore, u>(£) is negative for £ < 0, zero for £ = and is minimum for £ = £. For the stress tensor 
r{x,t) — —2u(x,t) we have (Fig. 3) : 



-d i r(0=m[p(0-p(0] K - w (0] 2 
i.e. the only points where = are £ = and £ = £. For the heat flux we have (Fig. 4) : 



mriQ (3 1 

~t~ w 



-K 2 



3 no 

2 T 



%q(0 = 2^)-p(0]K-M0] 



^K-^(C)] 5 



.«(0 

n(0 



and thus B^jq = 0for£ = 0,£ = £ and £ = £ where: 



r = - 



'3fc B r - 

m 



/ m 
KM 



If we adopt the thermal definition of temperature, i.e. fc^T = 2u/n, we find that (Fig. 5 a) : 

ddksTm = [p(£) - p(£)] {m[£ - <«(£)] 2 - k B T(0} 
and thus <9gT = if £ = £ and £ = £ where: 



(82) 

(83) 
(84) 

(85) 
(86) 



'feT - 



+ 



(87) 



Therefore d^T is negative for £ = a;/i g ]£,£[ and positive if £ ^ ]£,£[(Fig. 5 b). 

In conclusion with the thermal definition of temperature, the contribution jg9^(l/T) to the irre- 
versibility is non-negative if £ ^ ]£>£[ but negative for £ = x/t £ ]£, £[ . Therefore, whatever is the 
definition of entropy and of j$, it is impossible to obtain an equation of the form 



d t s + d x (sw+j s ) 



(88) 



with i — jQd^{l/T) without violating the Second Law (Fig. 6 a). Of course if we take Boltzmann's 
definition, Eq. (32), and js = Jb, then i(x, t) = 0. 
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7 Entropy and entropy production 

Whatever is the definition of entropy, from the Second Law 

d t s + d x {sw + j s ) = i, i(x,t)>0 (89) 
and the similarity of the flow, i.e. 



it follows that 
and 



s(x,t) = s(0, t.i{x,t)=t.i(£), £ = | (90) 
t.i(x, t) = t.i{£) = (w- £)<%s + sdjrw + d t j s (91) 

j t s~(t) = a de*.i(0-isK = o)} 

= I~ (t) — A js(C = 0) (92) 
If we take Boltzmann's definition of entropy, 

Sg{x,t) = 8g{() = 0[-£] |sg(-oo) + fc s y cfo [p~mp~ -p"lnp _ ] 

27rfc B T " 



1 



s b(-°°) = "0^2 



1 + ln 



m(n ) 5 



(93) 



then, as mentionned in section 3, if we define the entropy current by Eq. 31, i.e. js = js, we will 
have i(x, t) = and to obtain the Second Law in the usual form (with I~(t) = 0), i.e., 

I s - - WJ) P ^' )(t) (94> 

we have to introduce a "Boltzmann temperature" defined by 

T B (x,t) = (95) 

jB(X,t) 

However with such a definition we do not recover the usual phenomenological equations of thermody- 
namics. 

On the other hand if we take Boltzmann's definition of entropy and define the entropy current by 
the usual phenomenological equation jq = Tjs, where T = 2u/nks is the thermal temperature, then 
i(x,t) ^ given by Eq. (91). However from Eqs. (69) and (71) 

is negative. In other words, as already noticed in Section 6, i(x, t) can not be non-negative (see Fig. 6 
a) which is not compatible with the Second Law. 

If we adopt the definition of entropy given by the local equilibrium condition of a simple one- 
component fluid, i.e. 

.(0 = „o(o* fl i{i + m^!} (97) 



then 

As discussed in Section 6 this leads again to the conclusion that the Second Law is violated. 



T (0 = 7^§Y and /•:;-: j Q dJ±) ^ 
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8 Two- component fluid description 

To show that the time evolution of the piston problem satisfies the laws of thermodynamics we first 
remark that the fluid on the left side is characterized by 3 parameters (n^ , T ~ ,T P ). This leads us to 
look at our problem as a two-component fluid. 

At the surface x = 0, introducing 

ni(0,i) = ^, n 2 (0,t) = ^=, n = ni +n 2 (99) 

the Boltzmann entropy Eq. (67) can be express as 

s B (0,t) = s(u(0,t), ni (0,t),n 2 (0,t)) (100) 
s(u,m,n 2 ) = n k B { ^(1 +ln(4?r)) + ^ In - — r~ + g(n 1 ,n 2 )\ (101) 

and g(n 1 ,n 2 ) = C(5)=2 (— ^] In — + In ( 1 + — ^) - In 2 (102) 

\m+n 2 4 ) n 2 \ n 2 J 

We note that g(ni,n 2 ) — g(n\jn 2 ) implies 

ni-p- +n 2 -^ = (103) 
oni on 2 

From (101) we obtain for the temperature T s = du/ds and the pressure p, Eq. (8), at the surface of 
the piston: 

T -<»-'> = 2 » = T <°' f > = T ^ (1M) 

p(0,t) = n(0,t)k B T(0,t), i.e. r /r = r + p = 
where the subscript (+/— ) for right/left have been here omitted. 

In the bulk (x ^ 0) we shall assume that for all £ one can express the variables {no, To, Tp} in function 
of {u,ni,n2}, n = m + n 2 as was the case for x = 0. The thermal temperature will then be of the 
form 

T^) = ^-=T^,u,n u n 2 ) (105) 

and thus if we want to identify the temperature Eq. (105) with the thermodynamic temperature 
T s = du/ds the thermodynamic entropy must be given by 

= s{£,,u,n±,n 2 ) 

= n(0k B {l (l + ln^) +im-J+. g (£,n 1 ,n 2 )} (106) 

where at this point g{^,n\,n 2 ) is undetermined. For simplicity we choose g independent of £, thus 
given by Eq. (102), i.e. 

g{nun2 ) = 2 (^—-l\ ln^+lnfl + ^-ln2 (107) 
\ni + n 2 4/ n 2 V n 2 / 

Note that if n\{£) is chosen such that ni(£ = ±oo) = then we will have s = s B for £ = ±oo 
and for £ = ± 0, i.e. at those points where the fluids are in a stationary state. Moreover since 
g(ni,n 2 ) = g{n\/n 2 ) we have n\dgjdn\ + n 2 dg/dn 2 — which implies 

p = nk B T, i.e. r fr = T + p = 0. (108) 
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In that case it follows that 



. , | 1 \ , ■ r, ( fJ-2 ~ Ml 



fj, 2 - Mi = k B T 2 In h - — 

V ni 2\n 1 



which yields 



t.i(0 = -TjsjqdeT + 6k B ^ (0 £ n - 2^m) + o[gJ 
Moreover, the continuity equation for ni 

+9 € (ni«;+ii) = 

together with the boundary conditions 



lim ii(O = 0, 



lim m(0 = ^ 



implies (for £ < 0) 



ii(0 = - y) - n i w - 



Let us then look for a solution of the form 



(109) 



(110) 



(111) 



(112) 



(113) 



(114) 



(115) 



where at this stage Ai and A2 are unspecified constants. With this ansatz, and introducing 
x = — \//3£ > 0, we have for the left compartment 



7T 4/3 



/ m 



and 



£ (n - 2m) = 5n ^e- x2 {{x 2 - 1) [1 - An + A 2 (2a; 2 - 1)] - 2A 2 x 2 } + O 



/my 

\m) 



We thus have to find Ai and A 2 such that ti(£), Eq. (Ill), is non negative, where 



3Q(0&£ji = 62 "0** — (x 2 - l)(2x 2 - 1)[2 + x 2 (2x 2 - 1)] + O 



/my 



Moreover, we will have 



j t S~(t) = A £ d £ ti ft)-—^\j Q (p,t)+j 1 (to-»i)(-0,t)] 



(116) 
(117) 

(118) 
(119) 



with 



3Q(o,t) + ^(ja - m)(-o,t) = j H (-o,t) 



3qM 1--A 1V ^ 



In (l + 5") + 



1 5" 



2VTTT- 



3Q(0,t) 



1-5 —-^1^ 
16 



O 



(120) 
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and Jq(0) is given by Eq. (61). Using numerical computation, we see that choosing Ai = 1 and A2 = 0.1 
the entropy production Eq. (Ill) will be strictly positive (Fig. 6 b). For the right compartment (£ > 0) 
we have 



d_ 

It 



= £(m(0- y)-n 1 (0^)+^°°^(ni(0- y) (£ > 0) 

s+w = A y o ^ «co + r(+0; t) [jq (o, t) + ji ( M2 - mi ) (+0, t)} 



(121) 
(122) 



3Q(0,t) + Ml* - Mi)(+0,t) = jQ(0,t) 



1-5+^A^ 
16 



+ 



/ m \ 3 



(123) 



and Eq. (116-118) will remain valid with x = \J + £,- Therefore with Ai — 1 and A2 = 0.1 the entropy 
production i(x, t) is strictly positive for all (x,t), and all phenomenological equations of thermody- 
namics will hold (Fig. 6 a). 



Remarks 

1. With the choice Ai = 1 and A 2 = 0.1 the density m(x,t) corresponds approximately to the 
density n w (x,t) of particles with velocity larger than the average velocity of the fluid w(x, t) 
(see Fig. 7). 

2. We note that both definitions of heat current, (see Section 2 ) jq or jn — ]q + ji(^2 — Mi) 
arc quantities of order (m/M), which differ by a factor of order (m/M) 2 (see Eqs. (116-117), or 
(120)). However jq is continuous at the surface of the piston (as it should since it is the energy 
current), while jn is not continous (see Fig. 9). 

3. From Eqs. (116-117) follows that the possible choice of Ai, A 2 , leading to a positive entropy 
production, will be independant of the choice of the initial conditions (T^ , n J ) . In all the 
numerical computations (Figs. 1-9) the following parameters have been chosen 



1 



T - = l 



T, 



- 
i.e. 



10 



M = 10.000 
= 3000 



Po 



3000 



k B = 1 
300.085 
3000.85 



(124) 



4. Let us note that the difference between the entropy densities s(x,t) and sb(x, t) is very small 
of order (m/M) 2 (fig. 8 b). 



5. We recall that for one-dimensional ideal fluid the velocity of sound is y 3fc ^ T which coincide 
with Eq. (rcf85). 

6. If we consider the stationary states of our model, then the left compartment is an homogeneous 
fluid in a stationary state caracterised macroscopically by its temperature T (i.e. T~) its density 
n (i.e. n~), its velocity w(x, t) = and the heat flux jq. Microscopically this fluid is described 
by the density function 



p(x,v;t) 



0[- 



6[v] 2ni 



where, using the previous results, 

m n\ 



JQ 



2k B T n 2 



m n 2 
n 2 ni 



10 - 
— e 

7T 



l3v 2 



/ 

+ 6[-v] 2n 2 \ - e 

V 7T 







171 n 2 



2k B 
raw ' 



2k B T m 
n = n\+n 2 



(125) 



(126) 
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and the temperature of the piston is 

T P = T 



112 (M + m) 2 ( ri\ n 2 
rii AmM \ri2 n\ 



(127) 



The system is furthermore caracterised by its entopy function s = s(u, m, n 2 ), given by Eq. (106) 
and wc have: 

u = — n ksT p = n kpT = — r 

2 (128) 

k s[u(x,t),n 1 (x,t),n 2 (x,t)] = s B {x,t) 
Moreover in this stationary state ji = implies 



-S~(t) = -A J f 1 = (129) 

We obtain similar expressions for the right compartment and the stationary condition of the 
state is expressed by the relations 

P-=P + T P =T+ j q =j+ (130) 

Finally the total entropy production in the stationary state of this isolated system 

5* < 131 > 

is entirely due to the heat transfer from the subsystem at temperature T + to the subsystem at 
temperature T~which reflects the fact that the fluids are ideal. (Let us note that S = oo and 
thus only the entropy production is relevant in such non-equilibrium stationary states.) 

We should stress that the fluid in the stationary state is not locally in equilibrium (see Eq. (121)), 
which is a general property of far-from-cquilibrium stationary states, but is a two-component 
simple fluid. 



9 Conclusions and open problems 

Taking initial conditions such that on the average the piston remains at the same point so that there 
is only heat transfer from the right to the left compartment (T + > T^), but no work done, we have 
first derived the time evolution of the fluids on both sides of the piston. 

We have then shown that this fluid out of equilibrium can be considered as a 2-component simple 
fluid described by an entropy function s = s(u,n\,n2), and all the properties of the fluid will obey 
the laws of thermodynamics. In particular the two laws of thermodynamics are satisfied with 

- mnw 2 + u(s, m, n 2 ) 

du 2u 
ds nks 

sT + nifxi + H2^2 — u = nk B T = — r 

-Tw+j Q (132) 
where n, w, u, r, jg are defined by the usual relations of kinetic theory, /x, = du/drii 

Tjs=jQ+ji(l*2-Hi) (133) 



e = 

T = 

P = 

je = 



16 



and 



i = 3Q d x 



T 



+ jid x 



( M2 - Mi 



(134) 



is strictly positive. 

Moreover the entropy function which we have introduced coincides with the Boltzmann entropy 
in those domains where the fluids are in a stationary state. It was also shown that other natural 
definitions of entropy, temperature, irreversibility, which one could consider led to results violating 
the Second Law. Since the evolution is a similarity flow, one can not expect to obtain phenomenological 
relations of the form 



JQ 



C'u [j. 



C12 d x 



3i = C 2 i d x ( — ) + C 22 d x 



M2 - Mi 
T 

M2 - Mi 



(135) 



with Cij a state function (indeed d x 
However usual relations such as 



d 



dt 



s-(t) = i-(t) + 



,(+)-(-) 



T-(0,t) 



J+(t) 



P, 



(+)-(-) 



r+(o,t) 



with 
and 

arc satisfied with 



^(t) > 



K = A 



M 



1 

m7r 2 



> 



(136) 

(137) 
(138) 

(139) 



It is however an open problem to understand what is the physical meaning of the decomposition of the 
density as a sum of two densities. It is also to be understood whether there is a physical argument which 
would lead to a uniquely defined decomposition, and whether one could find a microscopical definition 
of our entropy function, i.e. whether it is possible to find a state function s(x, v; t) = s(p(x, v. t)) such 
that 



s(x,t) = s [u(x,t),ni(x,t),n 2 (x,t)] 



J dv s(x, v; t) 



(140) 
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Figure captions 



Figure 1: Density of particles n(x,t) — n(£) 
Figure 2: Fluid velocity w(x,t) — w(£) 
Figure 3: Pressure p(x, t) = p{£) = — r(£) 
Figure 4: Heat current j Q (x,t) = j Q (£) 

Figure 5: 

a) Temperature T(x,t) = T(£) 

b) Gradient <%T(£) 

Figure 6: Entropy production 

a) Contribution jg d^(l/T) 

b) i(x,t) = ti(£) 

Figure 7: Density ni(x,t)/n^ = rii(£)/"-o and 

density of particles n w (x,t)/n = n w (£)/n with velocity larger than w(x,t) 

Figure 8: 

a) Entropy density s ± (x, t)/s ± (±co) 

b) [s(x, t) — sb(x, t)]/s ± (±oo) 

Figure 9: Contribution [ji(n2 — (aj, t) = [ji(^2 — to the entropy current. 

Note that it is of the order 10~ 4 • ]q{x, t) 
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